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Abstract 

Goal, Scope and Background. The usefulness of power series 
expansion for an LCA system has often been doubted, as those 
systems may not possess the unique properties that enable power 
series expansion and analyses based on the power series. This 
paper surveys the existing literature on power series expansion 
of monetary input-output system and discusses how the power 
series expansion can be utilized for more general systems in¬ 
cluding the LCA model. 

Methods. The inherent properties of matrices that are capable 
of producing power series forms for their inverse and, further, 
can utilize structural path analysis are analyzed. Using these 
analyses, the way how a matrix that is not eligible for structural 
analyses is converted into an eligible form is investigated. A 
numerical example is presented to demonstrate the findings. 
Results. The necessary and sufficient condition for an indecom¬ 
posable, real square technology matrix can be expressed using 
power series was identified. Two additional conditions for a 
technology matrix to be utilized for structural analyses using 
power series expansion are discussed as well. It was also shown 
that an LCA system that fulfills the Hawkins-Simon condition 
can be easily converted into the form that is eligible for struc¬ 
tural analysis by rescaling the columns and rows. 

Discussion. As a numerical example, an application of accumu¬ 
lative structural path analysis for an LCA system is shown. The 
implications of the results are discussed in a more plain lan¬ 
guage as well. 

Conclusions. The survey presented in this paper provides not 
only the conditions under which a linear system is expressed 
using a power series form but also the way to appropriately 
convert a system to utilize the rich analytical tools using power 
series expansion for structural analyses. 

Recommendations and Perspectives. Widely used LCA databases 
and software tools have employed the linear systems approach 
as the basis. Much of these developments in the domain of LCA 
have been made, however, in isolation of the rich findings of 
IOA. There will be much to benefit LCA through an active dia¬ 
logue between the two disciplines. 

There are rich analytical tools available through the use of power 
series expansion. The current survey will help software devel¬ 
opers and LCA practitioners to apply such tools in LCA. 

Keywords: Input-output analysis; IOA; LCA practitioners; LCA 
software developers; monetary input-output system; power se¬ 
ries expansion; structural path analysis 


Introduction 

Doing an LCA study means working with numbers. The 
ecoinvent inventory data contains a few million numbers, 
and a typical impact assessment methods like Impact2002 
contains many thousands additional numbers. A bunch of 
numbers as such, however, makes no sense. They must be 
combined in the right way, using elementary operations like 
addition, multiplication, and division. The few sources that 
address the issue of combining the data in the right way are 
based on a formulation in terms of linear systems; see, e.g., 
Heijungs & Suh (2002) and Suh & Huppes (2005). 

As academic disciplines proliferate in the course of their 
developments, it is often the case that seemingly remote theo¬ 
ries in different disciplines in fact share much in common. 
Such recognition sometimes opens up the possibilities of 
gaining insights by using the findings from other disciplines 
or, sometimes, even of integrating them. The use of linear 
systems in economics, ecology and industrial ecology is an 
example that different disciplines share a common formal 
expression in representing their systems. Linear systems have 
been proposed, developed and used not only in LCA 
(Heijungs 1994), but in the fields of Input-Output Analysis 
(IOA; see Leontief 1936) and Ecological Network Analysis 
(ENA; see Hannon 1973) as well. But these developments 
have taken place in a rather isolated fashion. More recently 
there have been a number of efforts to utilize or merge the 
systems across the disciplines, using the common mathemati¬ 
cal ground as a basis (Joshi 1999, Suh 2004a, Suh 2004b, 
Suh 2005). The current paper concerns the linear system 
that is a generalization of those used in these fields its focus 
being on the conditions of power series expansion that en¬ 
ables various analyses. 

The use of power series expansion was introduced in the 
early history of IOA, where it was originally devised to solve 
the computational problem in inverting a large matrix un¬ 
der the lack of computational capacity (Waugh 1950). The 
power series form of a Leontief system has been extensively 
studied in relation to economic stability issues notably by 
Solow (1952). Although rapid development of computer 
technology quickly diminished the utility of power series 
expansion as a computational tool, its useful characteristics 
are still of analytical and computational interests (see e.g., 
Treloar 1997, Lave et al. 1995, Weber and Schnabl 1998, 
Lenzen 2001, Suh 2004a, Peters 2007). For instance, Treloar 
(1997) applied the power series expansion to extract impor¬ 
tant paths of the Australian inter-industry linkages in terms 
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of their embodied energy, and Lave et al. (1995) and Lenzen 
(2001) used it to quantify the truncation error of process- 
based LCA. Outside the Leontief input-output domain, 
Patten (1982) implicitly applied the power series form to 
analyze nutrient and energy flows between ecosystem com¬ 
ponents; Suh (2004a) used it to identify important suppliers 
in terms of their environmental impacts in a hybrid LCA 
setting. The virtue of using system expansion form is that it 
unravels a complex network relationship enabling a detailed 
insight on the structure of the system. 

Nevertheless, there have been confusions and misunderstand¬ 
ings on power series expansion especially for the systems 
other than the monetary Leontief input-output system, as 
those systems do not necessarily possess the unique proper¬ 
ties of Leontief input-output system that enables power se¬ 
ries expansion and analytical meaningfulness of the results. 
For instance, Frischkecht & Kolm (1995) and Schmidt (1995) 
applied the power series expansion form for LCA matrices, 
while their physical meaning is very doubtful (see Heijungs 
& Suh 2002, pp. 104-105). As it can be seen later in this 
paper, such a problem arises as physical systems, such as 
LCA and ENA models, may not share some of the key prop¬ 
erties of Leontief system that are related to power series 
form of an inverse. The same holds for Physical Input-Out¬ 
put Tables (PIOTs) and mixed-unit input-output tables where 
certain conditions are not met with these systems (Kratterl 
and Kratena 1990, Kratena et al. 1992, Pedersen 1999, 
Stahmer et al. 2003, Hoekstra 2003, Suh 2004b). Even within 
the monetary IOA, this situation may occur, for instance for 
heavily subsidized sectors (see footnote 1). 

This paper discusses how the power series expansion form 
can be utilized for a general system that includes LCA, mixed- 
unit input-output, PIOT and ecological network models. We 
discuss the conditions for a generalized linear system to have 
a power series form for its inverse, and, more importantly, 
the way to appropriately convert a system to utilize the rich 
analytical tools using power series expansion for structural 
analyses. Many LCA software tools and databases includ¬ 
ing some of the largest commercial packages have long been 
using linear systems for their embedded computational struc¬ 
tures. The linear systems used in these tools may not be eli¬ 
gible for structural analyses using power series form due to 
e.g. the tradition of treating waste flows. The current survey 
helps gaining insights on why some systems are not eligible 
for structural analyses and on how to manipulate such sys¬ 
tem to utilize the structural analyses. 

Many of the proofs or elements of proofs presented in this 
article are already made or implied by previous studies over 
the last century, e.g. Frobenius (1908), Oldenburger (1940), 
Waugh (1950), Solow (1952), Debreu and Herstein (1953) 
and Fiedler and Ptak (1962), but we felt it necessary to bring 
these rich findings from various other fields together, and, 
thus, they are included in this paper. The intention is to draw 
an attention to these findings, which are scattered over a 
long period of time, in a mathematically consistent way, 
rather than to develop original proofs. This then provides a 
basis for the development of new techniques in LCA, and 
for the incorporation of techniques that are already avail¬ 
able to specialists in IOA, ENA, or any other applied sci¬ 


ence that is based on linear systems, but that is not yet known 
to the LCA community. 

This paper is structured as follows. Section 1 discusses the 
power expansion for IOA, and Section 2 discusses more gen¬ 
eral systems, with an emphasis on LCA. Sections 3, 4 and 5 
show how the various results of Section 1 can be applied to 
LCA. An LCA application, using accumulative structural 
path analysis, in Section 6, demonstrates how techniques 
that were already available for IOA can now be employed 
to LCA systems. Section 7 concludes. At parts, the discus¬ 
sion will be rather technical, involving concepts from linear 
algebra like matrix norms, eigenvalues, and the Jordan ca¬ 
nonical form. Readers that are not familiar with such con¬ 
cepts are advised to move to the application, in Section 6. 
Or interested readers may find it useful to read a less techni¬ 
cal overview of linear systems used in LCA by e.g., Suh and 
Huppes (2005). 

1 Conditions for a Leontief System 

It was Waugh (1950) who, to our knowledge, first proved 
that a Leontief system always has an inverse that can be 
expressed as a power series. Let the direct requirement ma¬ 
trix be given as the square matrix B. The ;-/th element of a 
direct requirement matrix B, that is b tj , in an open Leontief 
system shows the amount of industry output from industry 
i directly required by industry j to produce one unit of out¬ 
put in monetary term. Waugh (1950) showed that the 
Leontief inverse (I - B) _1 can be expressed using a power 
series as 

(I - B) 1 = I + B + B 2 + B 3 + ... (1) 

provided that 

lim B 1 " =0 (2) 

Here we reproduce a proof following Chiang (1984, pp. 120- 
121). The equality in (1) is proved if 

(I - B) (I + B + B 2 + B 3 + B" ! ) = I, (3) 

where m—is assumed. The Left-Hand-Side (LHS) of 
(3) can be rewritten as 

(I + B + B 2 + B 3 + + B”’) - B (I + B + B 2 + B 3 + + B'") 

= (I + B + B 2 + B 3 + + B m ) - (B + B 2 + B 3 + + B'” +1 ), 

which becomes I - B'” +1 . Therefore, the condition that 
(3) holds if (2) holds. 

Waugh (1950) proved that a monetary Leontief system al¬ 
ways fulfills the condition shown in (2) using the norm of B. 
Let us define the norm of B as N(B), which here is the larg¬ 
est column sum of non-negative B (Bowker 1947, Waugh 
1950): ATB) = max(£.£..)_ 
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In general, a monetary Leontief system, representing a fea¬ 
sible economy, N(B) < 1, as, for each industry, the sum of all 
purchases from other industries to produce one unit of its 
output is generally less than unity given positive labor and 
profit. 1 One of the properties of a norm is that the product 
of the norms of two matrices is always larger than or at 
most equal to the norm of the product of the two matrices. 
For the same matrix, B, such property implies [N(B)] 2 
> N(B 2 ). As N(B) < 1 for non-negative B of a feasible 
economy, [N(B)] m always approaches 0 as m approaches 
positive infinity, which necessitates N(B m ) approaching 0. 
Hence, by the definition of a norm, all elements of B'” be¬ 
come 0 as m approaches positive infinity and thus the con¬ 
dition (2) always holds for any Leontief input-output sys¬ 
tem (Miller & Blair 1985). 2 

Solow (1952) provides a sufficient condition for a Leontief 
system to fulfill the condition in (2): all eigenvalues of a 
nonnegative, indecomposable Leontief system B must have 
a modulus less than unity. 3 There again, the condition that 
all column sums of B being less than unity is used to prove 
the theorem. 

It should be noted that the condition for all column sums of 
a Leontief direct requirement matrix being less than unity is a 
sufficient condition, not a necessary one. Furthermore, such a 
condition is not necessarily satisfied in the physical and mixed- 
unit systems such as the ones used for LCA, mixed-unit input- 
output tables and physical input-output tables (PIOT). 

2 Power Series Expansion for a more General System 

Linear systems are used for various fields outside the mon¬ 
etary input-output domain. For instance, LCA uses linear sys¬ 
tems to describe the commodity or functional flows between 
processes (see e.g., Heijungs 1994, Schmidt & Schorb 1995, 
Heijungs & Suh 2002, Suh 2004a), and many LCA databases 
and software tools including ecoinvent, Simapro 6 and 
CMLCA are employing such technique as a computational 
algorithm embedded in the database or the software tools. 

Mixed-unit input-output tables provide another example 
where linear system is employed. Konijn et al. (1997) and 
Hoekstra (2003) employed physical units for a number of 
metal flows in an input-output table, while other flows are 
noted in the traditional, monetary values. By using physical 
units, the input-output system can be free from the price 
inhomogeneity and fluctuation, better representing physical 
inter-relationships between industries. Another linear sys¬ 


1 1n reality it may be the case that some of the column sums of a technol¬ 
ogy matrix are larger than unity due to various reasons including that the 
capital investments are endogenized within the intermediate part, the 
sector actually made loss in that fiscal year or/and the sector is heavily 
relying on subsidies. 

2 We used the more direct proof by Miller and Blair (1985) instead of the 
one by Waugh (1950). 

3 Solow (1950) distinguishes decomposable and indecomposable matri¬ 

ces and provides a proof for each, and argues that the previous proof by 
Brauer (1947) is false, as it omits the decomposability condition. As the 
proof for a decomposable matrix is a directly corollary of the other, and 

as, in reality, it is hardly possible to have a decomposable matrix for the 
systems considered in this paper, we limit ourselves to indecomposable 
matrices throughout the paper, while the proof for the decomposable 
system is left to the readers. 


tem that deals with physical flows is the physical input-out¬ 
put table, PIOT. It shows the relationships between indus¬ 
tries in terms of the mass flows between them (Kratterl & 
Kratena 1990, Kratena et al. 1992, Pedersen 1999, Stahmer 
et al. 2003). Hubacek and Giljum (2003) utilize PIOT for 
analyzing the appropriation of land due to European ex¬ 
ports. Suh (2004b) borrowed the computational structure 
of LCA in dealing with the problem of representing waste 
flows in PIOT. 

Linear systems started to be used for ecosystems since 1970s 
(Hannon 1973). In ecology linear systems are used to de¬ 
scribe the structure of an ecosystem in terms of the energy 
and nutrient flows, and currently the world largest ecosys¬ 
tem network databases are based on the same, linear frame¬ 
work (Patten 1982, Szyrmer & Ulanowicz 1987, Christensen 
& Pauly 1992). Baily (2000) utilized the ecological network 
models for industrial systems and Suh (2005) generalized 
economic input-output model and ecological networks model 
in a unified framework. 

Although the system description in these applications looks 
very similar to the traditional Leontief input-output model, 
there are several notable differences. First, the most intui¬ 
tive difference is that, unlike the single unit monetary input- 
output tables, the column sum is not relevant in mixed-unit 
systems, as it adds, for instance, kg with MJ. Second, as a 
consequence, the norm of (I - A) is often bigger than 1, so 
that the sufficient condition for power series expansion in 
Miller and Blair (1985) and Solow (1952) is not necessarily 
satisfied in these systems. Third, some of these systems do 
not belong to the class of M-Matrices and some principal 
minors can be negative, while they still produce a non-nega¬ 
tive inverse (see Heijungs & Suh (2002), p. 17). 4 Fourth, 
LCA technology matrices are, at least as they are currently 
used, not necessarily used in a normalized form. Finally, the 
expression I - A has only physical meaning when the coeffi¬ 
cients on the diagonal of A are dimensionless, because 1 - a H 
contains for every i a dimensionless 1, and the algebra of 
units permits only the addition of quantities with the same 
unit. Thus these diagonal elements may be $/$, kg/kg, or 
MJ/MJ, but never $/kg, kg/MJ, etc. 

To illustrate these differences, consider, for instance, a simple, 
fictitious LCA system, consisting of a technology matrix A 
(2x2) and an intervention matrix E (1x2), as shown in Table l. 5 


Table 1 : Arbitrary LCA system, comprising 2 unit processes 



Power generation 

Coal mining 

Electricity (kWh) 

1500 

-3000 

Coal (Ton) 

- 0.01 

0.1 

CO 2 emission (kg) 

20 

10 


4 Fiedler and Ptak (1965) provide well-known equivalent definitions of M- 
Matrices and associated proofs, and Poole and Boullion (1974) offer a 
broader spectrum of M-Matrices with references. For example, a square 
matrix with non-positive off-diagonal elements with a non-negative in¬ 
verse is an M-matrix. 

5 Note that in this paper we followed the notation used in LCA, where the 
matrix A stands for what is generally noted as (I - A) in input-output 
analysis. To avoid any confusion, we used B for the direct requirement 
matrix of input-output accounts in this paper; see (1). Note also that use 
of A in this way is not common in input-output literature but is nothing 
new either (see e.g., Flawkins & Simon 1949). 


Int J LCA 12(6) 2007 


383 












Power Series Expansion 


Input-Output and Hybrid LCA 


Table 1 shows that producing 1500 kWh of electricity re¬ 
quires 0.01 ton of Coal, while coal mining produces 0.1 ton 
of coal by using 3000 kWh. An LCA technology matrix can 
be directly extracted from Table 1 such that 


and, thus, the matrix B'" has the eigenvalues /t/ w , ft 3 m , 
..., ft/”. Since both X -1 and X are independent of m, B m 
converges to 0 as m increases if and only if each of the p/" 
converges to 0. Therefore, B'” converges to a null matrix if 
and only if | ft | <1 for all i. 


1500 -3000 

"" - 0.01 0.1 ' 

The amount of total CO, emission is calculated by EA _1 = 
[0.0175 625], which shows that 1 kWh of electricity and 
lton of Coal generates 0.0175 kg and 625 kg of C0 2 emis¬ 
sion, respectively. This system works fine in producing its 
correct LCA result. 

However, although the system in (5) still generates correct, 
positive inverse 6 , its power series form, that is I + (I - A) + (I - 
A) 2 + (I - A) 3 + does not converge and thus all the useful 
analytical tools that make use of power series expansion 
form cannot be used for the system above. What follows are 
a discussion, on what conditions a more general system has 
a power series form for its inverse, and whether there is a 
way to convert a system that is not eligible for a power se¬ 
ries expansion like the one above to the one that has a power 
series form that can be used for various analyses. 


When B has less than n number of linearly independent eigen¬ 
vectors, X is not invertible and equation (7) does not hold. 
In that case, the Jordan canonical form can be utilized to 
complete the proof. For any B there exist a non-singular C 
that satisfies CBC -1 = J, where J is the Jordan canonical form. 
J has the following block diagonal form, 


J, o 
0 J 2 

0 0 


0 

0 


A 1 0 

0 ft 1 

0 0ft 

0 0 0 


0 

0 

0 


A 


A 1 0 

0 ft I 
0 0 ft 


0 

0 

0 


( 8 ) 


0 o 0 - ft 


3 Conditions of Power Series Expansion for a 
Generalized Linear System 

What follows is a rather technical elaboration on the condi¬ 
tions for a more general system to have a power series ex¬ 
pansion form for its inverse. These proofs will serve as a 
technical basis for subsequent discussions to follow. Non¬ 
technical readers may skip this section. 

Lemma 1: The inverse of an indecomposable, real square ma¬ 
trix A can be expressed as I + (I - A) + (I - A) 2 + ••• if and only if 
all eigenvalues of (I - A) have a modulus less than unity. 

Proof: Consider an nxn real matrix A. From (1) and (2), and 
by substituting A = I - B, it follows that the inverse of a 
matrix A can be expressed with a power series 

A- 1 = I + (I - A) + (I - A) 2 + (I - A) 3 + - (6) 

if and only if (I - A) to the power «z converges to a null matrix 
when m approaches infinity. An eigenvalue ft of B is defined as 
a possibly complex value that satisfies Bx, = ftx,. By locating 
eigenvectors x l5 x 2 , x 3 , ..., x„ to each column of X, BX = Xp 
holds, where is a square matrix in which ft, ft, ft, ..., ft, 
are on its diagonal and 0 on elsewhere. If X is non-singular, it 
follows that X~‘BX = p . By raising to the mth power 

(X -1 BX) m = (X _1 BX)(X~ ! BX)■ ■ ■ (X“‘BX) = X B X = p m 

- ■> (') 


6 The issue of positive inverse has been a point of discussion in input- 
output literatures especially in relation to the make and use framework 
and so called commodity assumptions (ten Raa et al. 1984, ten Raa 
1988, Kop Jansen & ten Raa 1990, Steenge 1990, Konijn 1994, Londero 
1999, Almon 2001). Negative elements in an inverse are better accepted 
in LCA and interpreted as the amount of inputs ’substituted 1 by non-pri¬ 
mary multiple outputs (see Heijungs & Suh 2002, pp. 45ff, Kagawa & 
Suh, forthcoming, cf. Stone et al. 1963). 


where /J,, ft,..., ft, are the eigenvalues of B, and all non¬ 
diagonal blocks are zero; the dimension of each Jordan block 
corresponds to the multiplicity of the eigenvector it repre¬ 
sents. By using the same technique used in (7), it follows 
that CB" ! C _1 “ = J”'. For m > ;, mth power of any jxj Jordan 
block, J/” is expressed by 


ft'" 

”'Qft" 1 

"C 2 ft- 2 - 

m c m . j+1 jur i+1 ~ 

0 

A," 

" Cft "- 1 ■■■ 

" H. A" ' 

O •• 

0 

A," 


• O 

0 

0 

ft" 


where the combinatorial factor n C r = «!/«! (n - r)\. By 
taking a limit, hm J, w iH vanish to a null matrix if and 
only if | ft | <1, and, thus, J"' will be 0 as if and only if 
| ft ] <1 for all i. 

Corollary 1: All real eigenvalues of a real, square matrix A 
that fulfils the condition in lemma 1 lies within the open 
interval (0, 2). 

Proof: Bx = (I - A)x = x - Ax = /ix <=> Ax = (l-/i)x = ko 
(1 - n) = X and, as | ft | <1 for all i if A satisfies the condi¬ 
tion in lemma 1, all real eigenvalues of A lies within the 
open interval (0, 2). 

Corollary 2: For a non-negative, real B the condition in 
lemma 1 is reduced to have the largest real eigenvalue of B 
within the open interval (0, 1). 

Proof: The Frobenius theorem adds a more powerful con¬ 
straint on the condition for a non-negative matrix, B. As the 
eigenvalue with largest modulus of non-negative B is real, 
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positive and simple (Frobenius 1908, p. 471), the necessary 
and sufficient condition for a non-negative, real B to be a null 
matrix by raising power of m , where m— »°° is that the larg¬ 
est real eigenvalue of B lies within the open interval (0, 1). 

Lemma 2: Let A be a square real invertible matrix, where 
B = I - A is non-negative. Then, the following statements 
are equivalent: 

a) (B) m = 0, as 

b) r c (B)<1, where r a represents the spectral radius, i.e. the 
modulus of the largest eigenvalue 

c) N(B m ) —> 0 

d) A -1 exists and is equal to I + B + B 2 + B 3 + ••• 

Moreover, a sufficient condition for any of these four state¬ 
ments is N(B)<1. 

Proof: See the previous proofs, Atkinson (1989, p. 490 
ff.) and Isaacson and Keller (1966, p. 14 ff. and 63). 

Lemmas 1 and 2 are well known and have been extensively 
discussed in input-ouptut literatures (cf. Hadamard's and 
McKenzie's theorems of dominant diagonal, see Takayama 
1974). Furthermore, these conditions hold for any arbitrary 
systems including LCA, PIOT, mixed-unit input-output sys¬ 
tem and ecological network models with a square matrix. How¬ 
ever, it should be noted that the lemmas present merely the 
condition for a system to have a proper power series expan¬ 
sion form for its inverse and does not tell us whether the power 
series form conveys analytical meaning to be used in, for in¬ 
stance, the Structural Path Analysis (SPA) and the Environ analy¬ 
sis. Such an issue will be dealt with in the following section. 


Now the question is whether the rescaled system in (10) can 
be used for structural analyses. Unfortunately applying A to 
the formula (6) gives meaningless figures instead of round- 
by-round requirements used in e.g., SPA. For instance, con¬ 
sider the third term, 


-0.5 

3 ' 

'-0.5 

3 ' 


'0.28 

1.2 ' 

0.01 

0.9 

0.01 

0.9 


0.004 

0.84 


This operation gives no information on upstream require¬ 
ments even though the power series form apparently con¬ 
verges to A -1 . 7 

In fact, the system requires one more condition in order to be 
utilized for round-by-round structural analysis: each column 
of A needs to be normalized by its diagonal element so that 
each column represents the input and output relationship per 
one unit of its output located in the diagonal. For instance, the 
system (10) will correctly give the upstream requirements when 
each column is normalized by its diagonal such that: 


1.5 -3' 

'1/1.5 

0 


1 

-30' 

-0.01 0.1 

0 

1/0.1 


-0.0067 

1 


Note that changing units was done by rescaling rows, while 
the operation in (12) rescales columns (see also Heijungs & 
Suh 2002. p. 109). Then, the second order of upstream re¬ 
quirement is calculated by 


4 Analytical Meaningfulness of Power Series Expansion 
for a more General System 

In a previous section it was shown that the system in (3) 
does not have a power series expansion form for its inverse, 
although the system correctly provides a solution. In fact, 
the eigenvalues of A are 1500 and 0.1, and therefore does 
not comply with the conditions of lemmas 1 and 2. 

By manipulating units, however, for instance kWh to MWh, 
(5) can be modified to have eigenvalues within the range 
specified in lemma 1 without changing the result. Consider 
a downscaled system of 


1.5 -3 

- 0.01 0.1 


( 10 ) 


0 

30' 

0 

30' 


'0.2 

0 ' 

0.0067 

0 

0.0067 

0 


0 

0.2 


(13) 


If interpreted, 1 MWh of power generation requires 0.0067 
ton of coal to produce one unit of its own commodity, and 
the 0.0067 ton of coal required by power generation, in turn, 
requires electricity, 30x0.0067 = 0.2 MWh. Therefore, in 
order to compute round by round requirements using a 
power series, the matrix should be manipulated in such a 
way that (net) output figure is located in its diagonal, and 
each column is normalised by its (gross) output figures. In 
case there is no self-consumption as is in many LCA sys¬ 
tems, this condition means 1 for all diagonal elements. 8 In 
addition, most applications will require non-positive off- 


The total C0 2 emissions calculated by using (9) are 17.5 kg 
and 625 kg per 1 MWh of electricity and 1 ton of coal, respec¬ 
tively, which are identical to the previous result. System (10) 
has norm bigger than 1 for (I - A), namely 3.9, which dis¬ 
qualifies the sufficient condition of lemma 2, but has eigenval¬ 
ues for (I - A) -0.521 and 0.921, and, therefore, fulfils the 
condition (b) of lemma 2. A general applicable type of rescaling 
is proposed by Janiszowksi (2003): the matrix A is replaced 
by aA with a chosen in such a way that the eigenvalues of aA 
are within the allowed range, e.g. a = 1/1 A max |, and the ob¬ 
tained inverse is multiplied by a. Thus, it is possible to use a 
power series expansion for finding the inverse of a general 
square matrix by introducing a suitable rescaling. 


7 For instance, using a plain flow diagram, the first order of upstream re¬ 
quirement on coal by electricity can be calculated by 0.01x3/0.1 = 0.3 
(ton coal/1.5 kWh electricity), while the RFiS of eq. (11) indicates 0.004 
(element in 2, 1). 

8 The self-consumption part appears in an aggregated statistics such as na¬ 
tional input-output tables. Unlike the common misunderstandings, this part 
does not represents the in-facility consumption of outputs such as the electric¬ 
ity used by the power plant that produces it, as what is generally counted is 
only the net output. These self-consumption parts in most cases show the 
inputs from other facilities that happen to be classified in the same industry 
category of the receiving facility. An example is food processing sector, where, 
for instance, a corn mill, classified as a food processing industry, produces 
starch and corn oil and these products are used by another food processing 
company for refinery and then by another to finally produce consumer prod¬ 
ucts. Such a statistical artifact is less likely to appear in LCA, where the pro¬ 
cess resolution is much higher. 
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diagonal, which assures that there will be only one com¬ 
modity produced by each process. 9 In other words, the sys¬ 
tem has to have been properly allocated 10 . Even if this addi¬ 
tional condition is not fulfilled, the round-by-round 
calculation can still be used, but its interpretation needs cau¬ 
tion. In this case, each round-by-round calculation will show 
the input requirements by that round minus the amount of 
input saved by producing the co-product that is assumed to 
be a by-product. In case the amount saved is larger than the 
amount used in that round, negative values will appear, which 
can be interpreted as the net amount of input saved by pro¬ 
ducing the co-product in that round. For a more detailed 
discussion that involves positive off-diagonal elements, see 
Lee (1982). 

The example shown above demonstrates that a system that 
fulfills the condition in Lemma 1 can be converted into the 
form that is eligible for structural analysis using power se¬ 
ries form for its inverse by rescaling its columns and rows. 
Then a natural question is whether such a rescaling of rows 
and columns alters the result or not. Let us consider a basic 
Life Cycle Inventory (LCI) problem, M = EA“*y, where E 
is environmental intervention per certain accounting period 
or operation time, A is a technology matrix based on the 
same accounting period or operation time, y contains the 
functional unit and M shows the inventory results (Heijungs 
& Suh 2002). Let us substitute A and y by fA and ry, re¬ 
spectively, so that all rows are arbitrarily rescaled. Then the 
RHS of the equation becomes EA _1 r _1 ry = EA”'y = M (see 
Heijungs & Suh 2002), so that the original solution is pre¬ 
served (for a seminal work, see Fisher 1965). Similarly, it can 
be shown that each column can be rescaled without altering 
the results. Let us substitute A and E by rescaled matrices, Ac 
and Ec> respectively. Then the RHS of the given equation 
becomes Ecc _1 A _1 y = EA _1 y = M, so that the result, M is al¬ 
ways preserved after the arbitrary rescaling of rows and col¬ 
umns. On top of that, it may be necessary to perform addi¬ 
tional operations on A, for instance changing the order of rows 
or columns. A general formulation is as follows: with trans¬ 
formations given by A' = RAC, E' = EC, and y' = Ry, the 
problem of finding a power expansion of A -1 thus reduces 
into the problem of determining matrices R and/or C, such 
that A'satisfies the conditions specified in Lemma 2. 

This series of discussion leads to the following : 

Corollary 3: The necessary and sufficient condition for a 
real, invertible A to be used for structural analyses using its 
power series form, A' 1 = I + (I - A) + (I - A) 2 + (I - A) 3 + 
is to satisfy all of the following three conditions: 

(a) a,-,- = 1, 

(b) a,y < 0 for all i, j ( i ^ /), and 

(c) all eigenvalues are within the unity in modulus for B = 
(I-A) 

The conditions (a) and (b) imply that the system is mono¬ 
functional, is arranged such that all production figures are 
located in its diagonal and that each column is related to the 

9 Having non-positive off-diagonal in a technology matrix implicitly applies 
so called 'by-product technology 1 model (see e.g., Konijn 1994). 

10 Allocation refers in LCA to the act of manipulating a technology matrix 

such that every process produces only one product. See Heijungs & 

Suh (2002) for more details and references. 


gross production value of the column. 11 Thus, a general 
recipe for constructing R and C can be made as follows: 

(1) first, appropriate steps must be made to transfer by-prod¬ 
ucts to other sectors, or to split sectors that produce more 
than one product in several sectors. Waste outputs should 
be changed into waste treatment inputs by changing the 
sign of an entire row of A and y. 

(2) we make R a permutation matrix, a matrix with a 1 in 
every row and column and zeros at all other places. The 
places where the 1 is to occur is determined by the dis¬ 
crepancy between the sector and 'its' product. 

(3) we make C a rescaling matrix that transforms the diago¬ 
nal of A into numbers between 0 and 1. 

In general, the make and use framework applied to a gen¬ 
eral linear system works this out (Kagawa and Suh, forth¬ 
coming). 12 

5 Relationship with Hawkins-Simon Condition 

In this section we discuss the relationship between the con¬ 
ditions for power series expansion of an inverse and the con¬ 
ditions for having positive inverse matrix. The Hawkins- 
Simon (H-S) condition states that a necessary and sufficient 
condition for a Leontief system B to have a positive inverse 
is that all principal minors 13 of (I - B) be positive (Hawkins 

6 Simon 1949). 

Let us generalize the H-S condition for an arbitrary system. 

Definition: If all principle minors of a real square matrix are 
positive, then the matrix is said to fulfill the generalized 
Hawkins-Simon condition. 

Note that, unlike the Leontief system, there exist some A 
that fulfils the generalized H-S condition, of which inverse, 
A -1 is non-positive. 14 However any A that fulfils the condi¬ 
tions (a) and (b) in Corollary 3 will always have positive 
solution (see, Fiedler and Ptak (1962)). Further, the general 
relationship between the condition in lemma 1 and the gen¬ 
eralized H-S condition can be easily derived. 

Corollary 4: Any square, real A that fulfils the condition in 
lemma 1 also fulfils the generalized Hawkins-Simon condition. 

Proof: The H-S condition is satisfied if and only if all real 
eigenvalues of A are positive, and, from the Corollary 1, all 
real eigenvalues of A that fulfils the condition in lemma 1 
lies within (0, 2). 15 

"In practice, column or row permutation matrices can be used to relocate 
the production figures to the diagonal. In case there should be a self¬ 
consumption portion due to e.g., poor statistical resolution, diagonals 
can be less than one but should still indicate the amount of net produc¬ 
tion per unit of gross production of itself (see also footnote 8). 

12 Models used under the make and use framework including commodity- 
technology, industry-technology and mixed technology models can be 
used for this purpose. However, for a mixed unit systems, certain re¬ 
striction applies: Industry-technology model requires a single unit for the 
part that the model is applied. Commodity-technology model does not 
impose such restriction. For an application of make and use framework 
for mixed-unit LCA system, see Suh and Huppes (2002). For a seminal 
review of make and use framework, see Konijn (1994). 

13 Forthe definition of a principal minor, we refer to Meyer and Mayer (2001). 
14 For example, the absolute form of A in (7), fulfills H-S condition, while its 
inverse is non-positive. 

15 The equivalence between the generalized Hawkins-Simon condition and 
all real eigenvalues being positive follows easily as the determinant of a 
matrix is the product of all eigenvalues of the matrix and the product of 
non-real eigenvalues of a real matrix is positive. 
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Lemma 3: For any real, square A that fulfils the conditions 
(a) and (b) of Corollary 3, the condition (c) is equivalent to 
the generalized H-S condition. 

Proof: As any A that fulfils the conditions in Corollary 3 will 
always have positive solution (see, Fiedler and Ptak (1962)), it 
is obvious that the generalized H-S condition is necessary con¬ 
dition for any A that fulfils (a) and (b) of Corollary 3 to be 
used for structural analyses. It remains to prove that Hawkins- 
Simon condition is a sufficient condition for any A that fulfils 
(a) and (b) of Corollary 3 to be used for structural analyses. 
For any A that fulfils (a) and (b) of Corollary 3 and a positive 
scalar c, a non-negative B exists such that A = d - B. Since, 
there exist a non-negative real eigenvalue /< max of non-negative 
B such that >\p\ for all eigenvalues p of B, and p = c - A, 
where A is the eigenvalue of A, ^ max is determined by A min as 
i^max = c where Vm is an eigenvalue of A such that A min 

< A for all A. Since A min is positive for all A that fulfils the 
generalized H-S condition 16 , | p | < c if A satisfies the general¬ 
ized H-S condition. Since a u < 1 for all A that satisfies (a) in 
Corollary 3, without any loss of generality one can set c = 1. 
Then, | p | < 1 for all eigenvalues p of B. Hence, by the defini¬ 
tion of generalized H-S condition, any A that fulfils the condi¬ 
tions (a) and (b) of Corollary 3 that fulfils the generalized H- 

5 condition also satisfies the condition (c) and vice versa. 

By the definition of generalized H-S condition and lemma 3, 
the lower bound of det(A) for A that satisfies the conditions of 
Corollary 3 is given. Furthermore, the upper bound of det(A) 
for A that fulfils the conditions of Corollary 3 can easily follow. 

Corollary 5: 0 < det(A) < 1 for any A that fulfils the con¬ 
ditions in Corollary 3. 

Proof: Since a a > 0 and < 0 (i^j) for all A that fulfils the 
conditions of Corollary 3, triangulizing A always adds non¬ 
positive values to each diagonal, whereas the result of 
triangulization will need to show all positive diagonal in 
order to have all positive principal minors. Therefore, each 
diagonal element of the triangulized A is equal or less than 

1. Since det(A) = (t 1 -t 2 -t 3 . t„), where t v t 2 , ... t n are the 

diagonal elements of the triangulized A, determinant of A 
that fulfils the conditions of Corollary 3 is less than unity. 

6 Application 

In this section, we demonstrate how power series expansion 
and corresponding analytical tools can be used for general¬ 
ized linear systems by means of a numerical example. The 
analytical tool that we chose is Structural Path Analysis (SPA) 
and Accumulative Structural Path Analysis (ASPA) based 
on Suh (2001). In analyzing a system it is often necessary to 
decompose the complex network system into individual 
paths. In this case the power series form of an inverse can be 
used to gain insights. Here we draw an example as an LCA 
problem, however, the methodology generally holds for other 
mixed-unit and physical flow models. Let us consider an 
LCA problem and define B = (I - A), then (6) becomes 

A -1 = I + B + B 2 + B 3 + - (14) 


16 See the proof of Corollary 4 and Fiedler and Ptak (1962) p. 385. 


For the given environmental matrix, E, the amount of envi¬ 
ronmental intervention per unit of output by each process is 
calculated by 

E (I + B + B 2 + B 3 + ■■■), (15) 

and especially the amount of environmental intervention i 
due to a unit output of process / is decomposed in scalar by 

e a + X e .A + X X e AA; +' ’' (16) 

l I m 

for all i, /, m, ... Note that each permutation of indices, /, 
m... in each term in (16) represents a specific input path 
required to deliver one unit of j to outside the system. A 
second-order of upstream path, XtAAs , for example, 
shows the amount of the zth environmental intervention gen¬ 
erated by the first commodity to produce fourth commodity 
to deliver one unit of the third commodity, which is set for 
the functional unit. 

With complete decomposition of an inverse shown in (16), 
SPA allows us to locate the key paths throughout the system 
that contribute significant environmental impacts. The most 
important paths from product p to process q with regard to 
the environmental intervention i is defined here as 

= (17) 

for all /, /, which means that the flow of product p to pro¬ 
cess q generates environmental intervention i the most among 
the single order paths. Similarly, the most important zzth order 
path, ( ■ "* r ) can be defined as 

n- \ n -1 

e i P A, V ’ •' Ax = .max (e u b tj b jm ■ ■ ■ b no ) > (18) 

for all indices /, 1, m , ..., o. The most important path in each 
order represents the most important packet of supply-chain 
at given length. In general, as the length of the supply-chain 
increases its impact decreases, since the 'contents 1 of far- 
upstream inputs in the final products are generally smaller 
than that of near-upstream inputs (Defourny & Thorbecke 
1984, Treloar 1997, Lenzen 2004) 17 . 

While SPA concerns only the direct generation of environmental 
intervention connected with each individual commodity path, 
ASPA takes both direct and indirect impacts into account. It 
identifies the bottleneck of a supply-chain where major envi¬ 
ronmental problems are caused by the particular commodity 
flow and the upstream supply-chain from that on. 

As an alternative to (16), the amount of environmental in¬ 
tervention i associated with a unit output of commodity / is 
decomposed in more compact form of scalar quantities by 

e i,+X e A’ ( 19 ) 

/ 

where a tj is z-/th element of A = A -1 . Or by 


17 A seminal review and analysis of available techniques can be found in 
Lenzen (2003). 


Inf J LCA 12(6) 2007 


387 







Power Series Expansion 


Input-Output and Hybrid LCA 


*» + X e A + XI>A.2w, (20) 

/ / m 

which is now expanded up to the third term. Likewise the 
equation (20) can be further expanded infinite times. Note 
that e l7 d (; in (19) shows the direct and indirect generation of 
environmental intervention i stemming from the first order 
input path of / to the ;th process including all upstream in¬ 
teractions originated from the particular input. 

In general SPA uses the certain rows and columns of the 
power series form in (6) and ASPA uses the power series 
form in (6) and the inverse, A -1 . For both analyses, it is ob¬ 
vious that the inverse needs to be expressed using power 
series form and, furthermore, the power series form needs 
to have the analytical meaning that we would like to utilize. 

Now we would like to use ASPA to analyze an LCA system. 
We utilized the ETH96 database that has been widely used 
for LCA studies (Frischknecht et al. 1996). ETH96 distin¬ 
guishes more than a thousand processes and several hun¬ 
dreds environmental interventions, focusing mainly on Eu¬ 
ropean product systems. It generally corresponds to the LCA 
computational structure by Heijungs (1994) and Heijungs 
& Suh (2002) that uses mixed-units. The technology matrix 
of ETH96 is transformed so that it fulfils the conditions of 
Corollary 3 as described in the previous sections, which in¬ 
volves assigning output values on the diagonals and nor¬ 


malizing each column by its diagonal. Beside, waste flows 
had to be modified as they are currently marked as outputs 
violating the condition (b) of Corollary 3. Therefore, the 
signs of all waste flows have been changed, so that the flow 
stands for 'waste treatment service' instead of the physical 
waste flows. The ASPA is carried out for the product, 'Elec¬ 
tricity, low voltage in UCPTE'. For the sake of simplicity we 
used C0 2 emission data as an indicator, while other indica¬ 
tors could have been be used as well. The result is shown in 
Table 2. Only those flows contribute more than 1% of the 
total are presented up to the third tier of upstream. 

Note that some of the commercial LCA software tools pro¬ 
vide a tree-like display, where the thickness of the arrows 
between the processes represents tier-wise accumulative con¬ 
tribution. In the existence of loops, i.e., when an output of a 
process requires itself through inputs from other processes, 
mathematical formulation of tier-wise accumulative contri¬ 
bution is not simple. For that, existing software tools with 
matrix approach utilize perturbation method to approxi¬ 
mate tier-wise accumulative contribution, which requires 
additional computational resources and time. 18 

The results show that 1.2% of the total C0 2 emission by 
the European low-voltage electricity mix are produced mainly 
for transmission within UCPTE, and almost all of the rest is 
made up by the electricity supplied by individual countries 
(1 st order). As to the contribution by individual countries, 


Table 2: Accumulative structural path analysis for ETH96 database (Electricity, low voltage in UCPTE) 



Order of upstream inputs 

% to the total 

I 8 ' 

2 n d 

3rd 

98.6% 

Electricity - Mix UCPTE 



1.2% 

Infra Electricity - in UCPTE 



42.0% 

Electricity - Mix UCPTE 

Electricity - Mix D 


20.2% 

Electricity - Mix UCPTE 

Electricity - Mix 1 


9.0% 

Electricity - Mix UCPTE 

Electricity - Mix E 


6.5% 

Electricity - Mix UCPTE 

Electricity - Mix NL 


6.4% 

Electricity - Mix UCPTE 

Electricity - Mix F 


5.2% 

Electricity - Mix UCPTE 

Electricity - Mix GR 


3.7% 

Electricity - Mix UCPTE 

Electricity - Mix B 


2.3% 

Electricity - Mix UCPTE 

Electricity - Mix P 


1.6% 

Electricity - Mix UCPTE 

Electricity - Mix A 


1.1% 

Electricity - Mix UCPTE 

Electricity - Mix Ex-YU 


1.0% 

Infra Electricity - in UCPTE 

Transmission 


20.4% 

Electricity - Mix UCPTE 

Electricity - Mix D 

Electricity from brown coal in D 

16.5% 

Electricity - Mix UCPTE 

Electricity - Mix D 

Electricity from hard coal in D 

10.6% 

Electricity - Mix UCPTE 

Electricity - Mix 1 

Electricity oil fired 1 

6.5% 

Electricity - Mix UCPTE 

Electricity - Mix 1 

Electricity gas fired 1 

5.8% 

Electricity - Mix UCPTE 

Electricity - Mix E 

Electricity from hard coal in E 

4.3% 

Electricity - Mix UCPTE 

Electricity - Mix GR 

Electricity from brown coal in GR 

3.4% 

Electricity - Mix UCPTE 

Electricity - Mix NL 

Electricity gas fired NL 

3.4% 

Electricity - Mix UCPTE 

Electricity - Mix D 

Electricity gas fired D 

3.3% 

Electricity - Mix UCPTE 

Electricity - Mix F 

Electricity from hard coal in F 

2.8% 

Electricity - Mix UCPTE 

Electricity - Mix NL 

Electricity from hard coal in NL 

2.7% 

Electricity - Mix UCPTE 

Electricity - Mix 1 

Electricity from hard coal in 1 

2.1% 

Electricity - Mix UCPTE 

Electricity - Mix B 

Electricity from hard coal in B 

1.9% 

Electricity - Mix UCPTE 

Electricity - Mix E 

Electricity from brown coal in E 

1.7% 

Electricity - Mix UCPTE 

Electricity - Mix F 

Electricity gas fired F 

1.3% 

Electricity - Mix UCPTE 

Electricity - Mix B 

Electricity gas fired B 

1.2% 

Electricity - Mix UCPTE 

Electricity - Mix P 

Electricity from hard coal in P 


UCPTE: Union for the Coordination of Production and Transmission of Electricity, W: Germany, I: Italy, E: Spain, NL: Netherlands, F: France, GR: Greece, 
B: Belgium, P: Poland, A: Austria, EX-YU: former Yugoslavia 
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electricity mix by Germany is the largest contributor account¬ 
ing for 42% of the total C0 2 emission by average low-volt- 
age electricity grid-mix in UCPTE, followed by that of Italy, 
contributing 20% (2 nd order). Regarding the German aver¬ 
age electricity, electricity from Brown coal (20%) and Hard 
coal (17%) combustion power plants dominates marking these 
two the most important third-order accumulated C0 2 em¬ 
bodiment paths. The largest CO,-contributing power tech¬ 
nology in Italy is oil-fired power plant, which is different from 
the case of Germany, accounting for 11% of the total (3 rd 
order). Using ASPA method the complex electricity grid-mix 
of Europe can be analyzed in terms of their contribution to 
certain emission. This example concerns only the C0 2 aspect, 
but of course, other indicators at different aggregation level 
such as total global warming potential, total human toxico¬ 
logical impact or total weighted impact can be used as well. A 
standard contribution analysis (Heijungs et al. 2005) of the 
same dataset does not specify the different orders, but effec¬ 
tively boils down to an infinite order ASPA. 

Likewise power series expansion and corresponding ana¬ 
lytical tools such as ASPA can be utilized in other physical 
systems including hybrid-unit input-output table, hybrid IO- 
LCA, PIOT and ecological network system, given that they 
fulfill the conditions discussed in this paper. 

7 Summary and Discussion 

In this paper, we discussed the conditions for an indecom¬ 
posable, real square matrix to have a power series form for its 
inverse in the context of structural analyses that are used in 
LCA, mixed-unit input-output analysis, PIOT and ecological 
network analysis. First, some known results for the monetary 
Leontief input-output system are reviewed. For the monetary 
input-output system, it is shown that the basic property of the 
direct requirement matrix B, namely the column sums of B 
are always less than unity, fulfills the sufficiency of having 
power series form for (I - B) _1 . For the physical systems such 
as LCA and mixed-unit input-output matrices, however, such 
sufficient condition is not met and, thus, the power series form 
may be unavailable. The condition for a physical system to 
have a power series form and, more importantly, the condi¬ 
tion for a system to be eligible for a structural analysis have 
been presented. The necessary and sufficient condition for an 
indecomposable, real square matrix A to have an inverse A -1 
that can be expressed as I + (I - A) + (I - A) 2 + •■• was identified 
to be that all eigenvalues of (I - A) have the modulus less than 
unity. In order to utilize structural analysis using power series 
expansion, however, A needs to fulfill two additional condi¬ 
tions, namely a„ = 1, and a l; < 0 for all i, j (i ^ /). It was also 
shown that a physical system that fulfills the H-S condition 
can be converted into the form that is eligible for structural 
analysis by rescaling the columns and rows. The relationships 
between different conditions are also discussed. Finally we 
demonstrated how a mixed-unit system can be transformed 
to utilize the analytical tools based on the power series expan¬ 
sion form using a numerical example. In the example we ap¬ 
plied ASPA for an LCA system to identify the most important 
accumulated path in the supply-chain. 


18 Suh (2005) presents a simple, scalar computation for the same issue in the 
context of Environ analysis used in systems ecology (Suh 2005, p. 262). 


In a more plain language, the condition in lemma 1 that all 
eigenvalues of a matrix (I - A) have a modulus less than 
unity simply requires (I - A) m approaches to a null matrix as 
m—>°° and vice versa. For this class of matrices, it is possible 
that the kind of physical interpretations used for round-by- 
round analysis is not given for the power series form of its 
inverse. The two additional requirements to be used for struc¬ 
tural analyses, that are a„ = 1, and a !; < 0 for all i, j (i ^ /), 
means that the system needs to be manipulated in such a 
way that each unit process, industry, compartment or eco¬ 
system component needs to produce only one flow, which 
needs to be located at the diagonal, and each column needs 
to be normalized by its gross production. In other words, 
each activity distinguished needs to be mono-functional and 
each column needs to present the inputs and outputs per 
unit of its gross output. In LCA such a manipulation can be 
done by allocation procedure, and for mixed-unit input-out- 
put model and PIOT, it can be done by using make and use 
framework. The generalized Hawkins-Simon condition im¬ 
plies that a system is not self-exhaustible. In other words, 
no activity or group of activities should require more than 
what is produced. The relationships between these condi¬ 
tions show that any mono-functional, normalized system 
(the conditions (a) and (b) of Corollary 3) that is not self- 
exhaustible (definition of generalized H-S condition) is eli¬ 
gible for structural analyses using power series form. 

The survey presented in this paper provides not only the con¬ 
ditions under which a linear system is expressed using a power 
series form but also the way to appropriately convert a system 
to utilize the rich analytical tools using power series expan¬ 
sion for structural analyses. Since Heijungs (1994), linear sys¬ 
tems have been widely utilized in LCA as a general database 
architecture for Life Cycle Inventory (LCI), as a basic computa¬ 
tional tool, and as a platform for wide range of analytical tools 
(see Heijungs and Suh 2002, for a seminal discussion). Widely 
used LCA databases and software tools such as Simapro 6, 
CMLCA, ecoinvent and its predecessor, ETH96 database all 
have employed the linear systems approach as the basis. Much 
of these developments in the domain of LCA have so far been 
made in isolation of the rich findings of IOA and ENA. As a 
result, the linear systems used in these tools may not be eli¬ 
gible for structural analyses using power series form due to 
e.g. the tradition of treating waste flows and the use of flow 
matrix instead of normalized coefficient matrix. The current 
survey tries to bridge the gaps between LCA on the one hand 
and IOA and ENA on the other hand, focusing on why some 
systems are not eligible for structural analyses and on how to 
manipulate such system to utilize the structural analyses. 

It is our conviction that a unified science, in which LCA, 
IOA and ENA are just different implementations of the same 
general systems-theoretical concept, is an ideal for which 
the LCA community should strive, together with other dis¬ 
ciplines of science. The subject of power series expansion 
elaborated here, however important, is only one example of 
how the rich findings of these other disciplines can strengthen 
the theoretical and analytical basis of LCA. 


Acknowledgement. The authors are thankful for the helpful comments 
by the two anonymous reviewers. 


Int J LCA 12(6) 2007 


389 








Power Series Expansion 


Input-Output and Hybrid LCA 


References 

Almon C (2001): Product-to-Product Tables via Product-Technology with 
No Negative Flows. Economic Systems Research 12, 27-43 
Atkinson KE (1989): An introduction to numerical analysis. John Wiley 
&C Sons, New York 

Bailey R (2000): Input-output modeling of material flows in industry. Doc¬ 
toral Dissertation, G.W. Woodruff School of Mechanical Engineering, 
Georgia Institute of Technology, Atlanta, GA, USA 
Bowker AN (1947): On the norm of a matrix. Annals of Mathematical 
Statistics 18, 285-288 

Brauer A (1946): Limits for the Characteristic Roots of a Matrix. Duke 
Mathematical Journal 13, 387-395 

Chiang AC (1984): Fundamental methods of mathematical economics. Third 
edition. McGraw-Hill Book Company, Auckland 
Christensen V, Pauly D (1992): ECOPATH II - a software for balancing 
steady-state ecosystem models and calculating network characteristics. 
Ecological Modeling, 169-185 

Debreu G, Herstein IN (1953): Nonnegative Square Matrices. Econometrica, 
21, 597-607 

Defourny J, Thorbecke E (1984): Structural path analysis and multiplier 
decomposition within a social accounting matrix framework. Economic 
Journal 94, 111-136 

Fisher FM (1965): Choice of Units, Column Sums, and Stability in Linear 
Dynamic Systems with Nonnegative Square Matrices. Econometrica 33, 
445-450 

Fiedler M, Ptak V (1962): On matrices with non-positive off-diagonal ele¬ 
ments and positive principal minors. Czechoslovak Mathematical Jour¬ 
nal, 12,382-400 

Frischknecht R, Bollens U, Bosshart S, Ciot M, Ciseri L, Doka G, Dones R, 
Gantner U, Hischier R, Martin A (1996): Okoinventare von Energiesys- 
temen. Grundlagen fiir den Okologischen Vergleich von Energiesystemen 
und den Einbezug von Energiesystemen in Okobilanzen fiir die Schweiz. 
Auflage No.3, Gruppe Energie-Stoffe-Umwelt (ESU), Eidgenossische 
Technische Hochschule Zurich und Sektion Ganzheitliche Systemanalysen, 
Paul Scherrer Institut, Villingen, www.energieforschung.ch, Bundesamt fiir 
Energie (Hrsg), Bern 

Frischknecht R, Kolm P (1995): Modellansatz und Algorithmus zur Berech- 
nung von Okobilanzen im Rahmen der Datenbank ECOINVENT. In: 
Schmidt M, Schorb A (eds), Stoffstromanalysen in Okobilanzen und Oko- 
Audits. Springer, Berlin, pp. 80-95 

Frobenius G (1908): Uber Matrizen aus positiven Elemente. Sitzungsberichte 
der Koniglichen Preussischen Akademie der Wissenschaften 26, 471-476 
Hannon B (1973): Structure of Ecosystems. Journal of Theoretical Biology 
41, 535-546 

Hawkins D, Simon HA (1949): Note: Some conditions of macroeconomic 
stability. Econometrica 17, 245-48 

Heijungs R (1994): A generic method for the identification of options for 
cleaner products. Ecological Economics 10, 69-81 
Heijungs R, Suh S (2002): The Computational Structure of Life Cycle As¬ 
sessment, Kluwer Academic Publisher, Dordrecht, The Netherlands 
Hoekstra R (2003): Structural Change of the Physical Economy - Decom¬ 
position Analysis of Physical and Hybrid Input-Output Tables, Doctoral 
dissertation, Free University of Amsterdam, Tinbergen Insitute Research 
Series, Amsterdam, The Netherlands 

Hubacek K, Giljum S (2003): Applying Physical Input-Output Analysis to 
Estimate Land Appropriation (Ecological Footprints) of International 
Trade Activities. Ecological Economics 44, 137-151 
Isaacson E, Keller HB (1966): Analysis of numerical methods. Dover Publi¬ 
cations, New York 

Janiszowksi KB (2003): Inversion of a square matrix in processors with 
limited calculation abilities. International Journal of Applied Mathematics 
and Computation Science 13, 199-204 
Joshi S (1999): Product Environmental Life-Cycle Assessment Using Input- 
Output Techniques. Journal of Industrial Ecology 3, 95-120 
Kagawa S, Suh S (forthcoming): Multistage process-based make-use system. 
In: Suh S (ed), forthcoming: Handbook of Input-Output Analysis for 
Industrial Ecology, Springer 

Konijn PJA (1994): The make and use of commodities by industries-On the 
compilation of input-output data from the national accounts. Doctoral 
dissertation, University of Twente, Twente, The Netherlands 
Konijn PJA, de Boer S, van Dalen J (1997): Input-Output Analysis of Mate¬ 
rial Flows with Application to Iron, Steel and Zinc. Structural Change 
and Economic Dynamics 8, 31-45 

Kop Jansen P, ten Raa T (1990): The choice of model in the construction of input- 
output coefficients matrices. International Economic Review 31,213-27 
Kratena K, Chovanec A, Konechy R (1992): Eine okologische volkswirt- 
schaftliche Gesamtrechnung fiir Osterreich. Die Umwelt Input Output 


Tabelle 1983. Institute fiir sozial-, wirtschafts- und umweltpolitische For- 
schung, Wien, Austria 

Kratterl A, Kratena K (1990): Reale Input-Output Tabelle und okologischer 
Kreislauf. Physica-Verlag, Heidelberg, Germany 
Lave L, Cobas-Flores E, Hendrickson C, McMichael F (1995): Using input- 
output analysis to estimate economy wide discharges. Environmental 
Science St Technology 29, 420-426 

Lee KS (1982): A Generalized Input-Output Model of an Economy with En¬ 
vironmental Protection, Review of Economics and Statistics 64, 466-73 
Lenzen M (2001): Errors in conventional and input-output-based life-cycle 
inventories. Journal of Industrial Ecology 4, 127-148 
Lenzen M (2003): Environmentally important paths, linkages and key sec¬ 
tors in the Australian economy. Structural Change and Economic Dy¬ 
namics 14,1-34 

Leontief WW (1936): Quantitative Input and Output Relations in the Eco¬ 
nomic Systems of the United States, The Review of Economic Statistics 
18,105-125 

Londero E (1999): Secondary products, by-products and the commodity 
technology assumption. Economic Systems Research 11, 195-203 
Meyer CD, Meyer C (2001): Matrix Analysis and Applied Linear Algebra, 
SIAM, Philadelphia 

Miller R, Blair P (1985): Input-Output Analysis: Foundations and Exten¬ 
sions. Prentice-Hall, Englewood Cliffs, NJ, USA 
Oldenburger R (1940): Infinite Powers of Matrices and Characteristics roots. 

Duke Mathematical Journal 6, 357-361 
Patten BC (1982): Environs - Relativistic Elementary-Particles for Ecology. 
American Naturalist 119, 179-219 

Pedersen OG (1999): Physical Input-Output Tables for Denmark. Products 
and Materials 1990, Air Emissions 1990-92. Statistics Denmark, Copen¬ 
hagen, Denmark 

Peters G (2007): Efficient algorithms for life cycle assessment, input-output 
analysis, and Monte-Carlo analysis. Int J LCA 12 (6) 375-382 
Poole G, Boullion T (1974): A Survey on M-Matrices. SIAM Review 16, 
419-427 

ten Raa T, Chakraborty D, Small, JA (1984): An alternative treatment of 
secondary products in input-output analysis. Review of Economics and 
Statistics 66, 88-97 

ten Raa T (1988): An alternative treatment of secondary products in input-output 
analysis: Frustration. Review of Economics and Statistics 70, 535-538 
Schmidt M, Schorb A (eds) (1995): Stoffstromanalysen. In: Okobilanzen 
und Oko-Audits. Springer, Germany 

Solow R (1952): On the structure of linear models. Econometrica 20,29-46 
Stahmer C, Kuhn M, Braun N (2003): Physische Input-Output Tabellen, 
Beitrage zu den Umweltokonomischen Gesamtrechnungen. Metzler- 
Poeschel Verlag, Stuttgart, Germany 

Steenge AE (1990): The commodity technology revisited-Theoretical basis 
and an application to error location in the make-use framework. Eco¬ 
nomic Modeling 7, 376-387 

Stone R, Bacharach M, Bates J (1963): Input-Output Relationships, 1951- 
1966, Program for Growth. Vol. 3. London, Chapman and Hall, UK 
Suh S (2001): Accumulative Structural Path Analysis for the U.S., Presented 
at the Danish Prioritization Study Workshop, Copenhagen, Denmark 
Suh S (2004a): Functions, commodities and environmental impacts in an 
ecological economic model. Ecological Economics 48, 451 - 467 
Suh S (2004b): A Note on the Calculus for Physical Input-Output Analysis 
and its Application to Land Appropriation of International Trade Ac¬ 
tivities. Ecological Economics 48, 9-17 
Suh S (2005): Theory of Materials and Energy Flow Analysis in Ecology and 
Economics. Ecological Modeling 189, 251-269 
Suh S, Huppes G (2002): Economic input-output models for allocation in 
LCA, SETAC-Europe annual meeting, Wien, Austria 
Suh S, Huppes G (2005) Methods in Life Cycle Inventory (LCI) of a prod¬ 
uct. Journal of Cleaner Production, 13 (7) 687 - 697 
Szyrmer J, Ulanowicz RE (1987): Total Flows in Ecosystems. Ecological 
Modeling 35, 123-136 

Takayama A (1985): Mathematical Economics. 2 nd edition, Cambridge, UK, 
Cambridge University Press 

Treloar G (1997): Extracting embodied energy paths from input-output 
tables: towards an input-output-based hybrid energy analysis method. 
Economic Systems Research 9, 375-391 
Waugh FV (1950): Inversion of the Leontief Matrix by Power Series. Eco¬ 
nometrica 18,142-154 

Weber C, Schnabl H (1998): Environmentally Important Intersectoral Flows: 
Insights from Main Contributions Identification and Minimal Flow 
Analysis. Economic Systems Research 10, 337-56 

Received: January 27th, 2006 
Accepted: August 23rd, 2007 
OnlineFirst: August 24th, 2007 


390 


Int J LCA 12 (6) 2007 



